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1.0 INTRODUCTION 


The concept of the "low- low” GRAVSAT satellite to satellite tracking mission 
has been proposed as a technique for helping to determine the fine structure of 
the earth’s gravity field. The system consists of txro drag-free satellices 
separated by a few hundred kilometers in the same low altitude polar orbit con- 
strained to follow a purely gravitational path. Such a satellite configuration 
is highly sensitive to geopotential variations, and the polar circular orbit 
would provide continuous tracking over the entire earth, with either one or both 
satellites having tracking capability and ability to transmit telemetry to the 
earth. The tracked satellite would have an ordinary spacecraft transponder. 

The physical reason for sensitivity to local gravity perturbations in the 
"low- low" configuration is clear. As the two satellites approach a gravitational 
perturbation, the nearer satellite "feels" its presence sooner and experiences an 
acceleration relative to its sister satellite. When the first satellite passes 
the perturbation, its acceleration changes sign relative to its sister, and as it 
recedes from the perturbation, the acceleration it "feels" diminishes. The sister 
satellite experiences a similar acceleration, but out of phase with the first 
satellite. The projection of the difference between their velocities along the 
separation vector, the satellite-to-satellite range rate, then has a definite 
signature for the satellite configuration and the gravity pertiirbation . Gravity 
perturbations far from the two satellites tend to act on each satellite in the. 
same manner and hence produce little relative velocity (Figure 1) - 

Highly accurate gravity models are needed for two purposes ' to accurately 
predict satellite motion and to determine the detailed shape of the geoid. NASA’s 
Earth and Ocean Physics Program (EOPAP) will require a deter m ination of the geoid 
to 10 cm, together xrlth an orbital prediction accuracy in satellite altitude to 
10 cm. The requirement placed on the accuracy of the gravitational fine structure 
necessitates an extremely large number of parameters in the gravity field model. 

A simultaneous data reduction for all model parameters would be a tremendously 
difficult computer task, even with the most sophisticated new computer systems. 

The range rate data between the two satellites in the "low-low" configuration, 
however, is highly sensitive to local gravity perturbations and much less sensi- 
tive to gravity perturbations distant from the satellite pair. This supports the 
attractive idea of representing the gravitational fine structure in a local 


1 



Satellite #1 



time -> 




The top two curves diagram the change in velocity magnitude 
as a function of time as the low-low satellites pass over a 
perturbing anomaly. The bottom curve is the difference of 


the top curves , showing the approximate shape of the relative 
range-rate signal between the spacecraft .caused by the 
perturbation- 


FIGURE 1 



mathematical model and estimating a local area of gravitational perturbations 
using local data. This would, of course, require a fewer number of model para- 
meters and make the 'problem more manageable. The global solution could then be 
built up from a set of these local solutions- This idea has been tested by 
simulation for SST data by Schwarz [9] and Hajela [4]. 

Schwarz utilized specially developed software for least squares parameter 
estimation to study the feasibility of the "low-low"- and "high- low" mission with 
summed SST range and range-rate data, i.e. the measurement proceeded from a ground 
station to the tracking satellite to the tracked satellite and back again. He 
selected surface density blocks to represent the high order gravity field above 
a truncated spherical harmonic reference field and limited his study to a geo- 
graphically isolated set of these blocks. Data was simtilated over the blocks and 
a least squares estimating process attempted to recover the surface density para- 
meters. His study .processed only short data arcs and solved for only the surface 
density block parameters. The orbit states were considered to be nuisance para- 
meters and were not estimated. A serious shortcoming of his study was neglecting 
the contaminating effects of errors iii non-estimated neighboring surface density 
blocks on the accuracy of the’ estimated blocks. We refer to this as aliasing. 

Hajela, using the GEODYN program [2] performed a similar feasibility study 
for the "high-low" mission using equal area' gra-vity anomaly blocks to represent 
the high order gravity field. Again using simulated swmted SST range and range- ' 
rate data he attempted a least squares recovery of local sets of gra-vity anamoly 
blocks using multiple short arcs of data. Unlike Schwarz, he included the alias- 
ing effects of errors in neighboring blocks on the estimated blocks by Including 
in his estimated solution a band of blocks around the perimeter of the blocks of 
interest and outside of the data coverage. In his study, Hajela essentially 
assumed that he knew the satellite states perfectly by assigning them very small 
a priori sigmas in the estimation algorithm. T^hereas Schwarz performed all of 
his least squares adjustments using infinite a priori sigmas for the surface 
density blocks, Hajela included in his investigation cases where a priori informa- 
tion for the gravity anomaly blocks was assumed. 

The purpose of this investigation is to apply covariance error analysis 
techniques to determine under what estimation strategies the "low-low" mission is 


3 



feasible for local xecovexy of gravitational fine structure when the effects of 
aliasing are considered. This investigation treats 5° x 5° surface blocks with 
the satellite altitudes at 250 km. Other studies (Schwars [9], Lowrey [8]) have 
investigated the sensitivity to and resolvability of surface features of given 
size as a function of the "low-low" orbital altitude. ‘ The drag compensation 
system of the GRA.VSAT satellite has been assumed to be error free in this study, 
so that the results obtained xjill be somewhat optimistic. In a more detailed 
and realistic study the aliasing effect in the parameter estimation due to imper- 
fect knowledge of surf ace . force compensation must be included. 


<®’ POOR 
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2.0 SUR1?ACE DENSITY GEOPOTENTIAL REPRESENTATION 

The purpose for the present study of the low-low SST experiment is to deter- 
mine the feasibility of recovering the fine structure of the earth's gravity 
field as a composite of recoveries over local regions. The total geopotential is 
represented as 

W = U + T 

where the function U is assumed to be a fixed low order reference field and T 
represents the high order potential. 

The potential T could be expressed by any one of several mathematical models, 
i.e. spherical harmonics, gravity anomalies, saiE5>ling fmctions, spherical grid, 
etc. Whatever model is chosen, the description of the gravity fine structure 
will require a very large set of parameters. For example, 2 ° geoid resolution 
would require in excess of 8400 spherical harmonic coefficients. While the "low- 
.low" SST experiment will provide global data coverage and the simultaneous . esti- 
mation of the full parameter set- is theoretically possible, it represents a very 
large and difficult numerical task for even the largest computers, ^^hile a 
complete study of such a global fine structure estimation is worthwhile, it will 
not be pursued in this investigation. Instead we attempt to estimate subsets of 
parameters while keeping the remainder at a priori values. The complete solution 
is then built up from these subset solutions. However, mless the parameterization 
exhibits local independence or "orthogonality" with respect to the data type, the 
uncertainties in the neglected parameters will badly alias the adjusted parameters 
and such a recovery procedure will not be possible. Orthogonality has been defined 
by Argentiero [1] as follows: For a given estimation procedure, the jth adjusted 

parameter is orthogonal to the kth unadjusted parameter if the aliasing contribu- 
tion of k to j is zero. 

We have selected to parameterize the high order geopotential by local surface 
density blocks representing a fictitious surface density layer. 
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where G is the gravitational constant, da is the element of surface area, p is 
the density of the mass layer and d is the distance from the point where T is 
evaluated to the integration element da. The parameter used to describe the 
high order field is ® = Gp with units of acceleration. The value of the parameter 
$ will be expressed in milligals (mgal) , where 

—3 —3 2 

1 mgal = 10 gals = 10 cm/sec 

The relationship of the parameter $ to gravity anomalies Ag is given approximately 
by (Heiskanen and Moritz, page 303) 

$ = -i- (Ag + — g — ) • 

2tt ^ ^ 2 ^ . 

where g .is the mean normal gravity, N is the geoid undulation and R is the mean 
earth radius. For order of magnitude comparisons' between the surface density 
parameter ^ and corresponding gravity anomalies, the expression 

$ = As 

2tT • - - ■ ; 

is adequate. A detailed treatment of the surface layer potential is given by 
Schwarz, including transformations between the density layer representation. and 
spherical harmonic and gravity anomaly representations. For simplicity, the mass 
layer is considered spread on the surface of a sphere with radius equal to the 
mean radius of the earth. - • , ' 

The density layer is modeled by individual area blocks in which the density 
is a constant value representing the average density of the surface layer within 
the block. The disturbing potential is then 

where the integration is over the area of the jth block and the sum is over all 
surface blocks. The integral over each block is evaluated by dividing the block 
into a number of sub-blocks of area a,, , where the distance d , is the distance 

jk jk 

from the center of the sub-block to the point of evaluation, ' . ' 

OEIGINAI] PAGE IS. ' 
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Note that this is equivalent to representing the potential of each block by the 
sum of the potentials of a point mass located at the center of each sub-block, 
with mass The accuracy with which this summation of sub-blocks represents 

the integral over the block depends, of course, on the number of sub-blocks used. 
For satellite distances which are large compared to the block size, a single mass 
point approximation is suitable, whereas for satellite disturbances which are of : 
the same order as the block size, many mass points are needed (Figure 2) . 


The principal portion of this investigation is concerned with recovering local 
subsets of blocks utilizing surf ace. blocks with boundaries defined by constant 
values of latitude and longitude increments (equalangular blocks). Such blocks 
have areas which are strongly dependent on latitude, e.g. a 5° x 5° block on the 
equator possesses a larger area, and hence produces a larger perturbation for a' 
given surface density, than a 5° x 5° block at a different latitude. 
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3.0 COMPUTER SOFTWARE DESIGN 


In. an effort to utilize existing software in fhe investigation, where possible, 
several existing orbital estimation /error analysis computer programs were con- 
sidered: GEOUra/ERODYM, NAP/NAPCOV, ORAN and GTDS. Constraints placed on the 

ultimate choice of software were: 

1) Adequate documentation. • 

2) Capability for a large set of adjusted parameters. 

3) Multi- arc/multi-satellite processing. 

4) SST data type. 

5) Surface density force model. 

6) Data generation capability. 

7) Parameter estimation capability. 

8) Error analysis capability. 

9) Availability of program for modification. 

Although none of the systems available -filled all needs, the program most 
appropriate was determined to be GEODYN/ERODYN. 

GEODYN is a large scale orbital parameter estimation program oriented toward 
geodynamlcal applications, while ERODYN is an orbital and geodetic error analysis . 
program designed to be operated directly from GEODYN output. To perform an error 
analysis of the GEODYN estimation algorithm, ERODYN requires the estimation 
normal matrix output from GEODYN. Moreover, to investigate by error analysis the 
partitioning of a parameter set into adjust and unadjust (alias) parameters, the 
total parameter set must be adjusted in a GEODYN run so that all parameters are 
included in the normal matrix. This normal matrix, obtained from a single GEODYN 
run, is then partitioned into appropriate adjust and unadjust sets within ERODYN. 
The setting of a priori parameter sigmas and the partitioning of the total para- 
meter set into adjust and unadjust parameters within ERODYN permits a great many 
error analysis runs to be made from a single GEODYN normal matrix. 

As the surface density gravity model and the SST data type capability were 
a little used portion of the GEODYN pro gram^ extensive checking of these features 
was performed and program errors corrected. The variational equations were 
checked by numerical secant partials, i.e. by incrementing the initial state and 



surface density parameters and integrating the state vector so’;that partial 
derivatives could he formed relative to the nominal state vector. 

The relative range and range-rate measurements hettvreen the ttw satellites and 
the measurement partials are calculated by the program in terms of the difference 
between the state vectors of the satellites and the geometric and variational 
partials 



where the state vectors and variational partials are integrated separately for 
each satellite. Calculating relative motion quantities in such a manner can lead 
to important significant digits being lost in the subtraction. As a check on the 
ntimerical accuracy maintained, a closure test was performed on both the satellite 
state vector and variational partials. By closure we mean the integrating of the 
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system of differential equations forward in time from an initial solution at 

time t to a solution X, at time t^, and then reversing the time direction and 
o It 

integrating backward from the solution X^ at time t^ to a solution X^ at time 
t^. A comparison of X^ and X^ reveals the accuracy of the integration process. 

Our simulation for satellite orbits of interest for the "low-low" mission using . 
an 11th order Cowell integrator disclosed that 12 significant digits in the 
satellite state and 9 digits in the variational equations were retained when 
t^ “ t^ was set at 8 hours and an integration step size of 20 seconds was used. 

This accuracy assures 4 to 5 significant digits in calculated relative motion 
quantities, which is adequate for error analysis purposes. However, for future 
estimat'ion/TecoveT^ purposes using relative range or range-rate SST data, it will 
probably he necessary to integrate the system of orbital equations for the relative 
state vector, where appropriate rearrangement has been performed to eliminate the 
subtraction of nearly equal quantities [3],’ [7]. 

GEODYN, as indicated, is primarily a scientific tool for geodetic and geo- 

dynamical study and is not easily implemented for mission analysis needs. In 

particular, there is no efficient data simulation capability in the system. To 

fulfill the needs of this study significant modifications to the GEODYR system 

■( 

version 7410.0 were performed, rendering it a more applicable mission study soft- 
ware system. This version of GEODYN with the mission analysis capability updates 
has been denoted GEOllAP, Geo d ynamics and Mission ^alysis 'program. Some of the 
modifications to the GEODYN/ERODYM system for the performance of this study 
include : 

? SST range, range- rate, summed range and summed range-rate data generation. 

e Selection of simulated data by data t 5 npe, time intervals, and geographic 

latitude and longitude limits for the ground track. 

• Equal- area surface density blocks. 

«• Additive random noise for simulated data. 

“ Simplified user input for surface density blocks. 

• Simplified force model evaluation for surface density blocks. 

• Multi- arc capability for ERODYN. 

® Triple precision matrix inversion routine for ERODYN. 
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• HCL (radial, cross track, along track) state vector a priori sigma 
input for ERODYN- 

Modification to the force model evaluation for surface density blocks was made 
in an effort to reduce computer time. The effect of a surface density block is 
modeled in the program by a user specified number of mass points symmetricaliy 
placed within the block. However, only blocks near the satellites need such 
accurate mathematical representation and blocks far from the ground track may be 
modeled much more simply. GEOMAP has been implemented to offer three options: 

1. All surface blocks represented by the same number of user specified 
mass points. 

2. All surface blocks within a cone defined by a user specified geocentric 
angle from the sub satellite point represented by the user prescribed 
n nmb pi-r of mass points , and blocks outside of this angle, represented by 
a single mass point. 

3. All surface blocks outside a given angular cone of the subsatellite 
block completely neglected. This option may be used in conjxmction 
with options 1. or 2. 

1 

Option 2 is similar to the technique implemented in later versions of GEODYN 
where the number of mass points representing the blocks inside of the cone is set 
to 9 and is not a user input. 

As implemented at the time of this study, ERODYN operated as a single arc 
error analysis package. Multi- arc capability was added by modifying a program 
made available by Dr. P. Argentiero (private commtinication) designed to concatenate 
single arc normal matrices into a single larger dimensional normal matrix suitable 
for analysis by the linear algebra capability of ERODYN. This leads to some 
degree of inefficiency in that blocks of zeros are introduced into the concatenated 
normal matrix which are subject to mathematical computations. However, as a means 
to quickly add multi-arc error analysis capability to ERODYN for the purposes of 
this study it is entirely adequate. 

In light of the modifications made to GEODYN Version 7410.0 and to avoid 
possible confusion, the software system used for this study will be designated 
GEOMAP /ERODYN in the remainder of this report. 
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4.0 COVARIANCE ERROR ANALYSIS DESCRIPTION 


The error analysis for the Bayesian least squares estimator in GEODYN is 
described in the ERODYN Program Mathematical Description [5]. However, for 
completeness xre reproduce the main points of the mathematical treatment. 

Let y be an n dimensional vector of observations modeled by the nonlinear 
equation 


y = E(X,S) + e 


V7here X and S are parameters of the dynamical system. Here X denotes the para- 
meter set to be adjusted (estimated) and S denotes the parameter set which is 
unadjusted, or constrained to constant values in the estimation process but whose 
uncertainties are to be considered in computing the covariance matrix of the 
estimate. The vector e represents the observation noise, which is mean zero and 
statistically Independent. Linearizing about a nominal solution of the parameters 
2^ and yields 


r<5Xl 


fiy = C 


+ e 


SS 


= A6X + B(SS + £ 


where 


5y = y - P(x^,s^) 

6X = X - Xj^ 

6S = S - S„ 

N 


A = 


8F(X,S) 

8X 
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and 


B = 


3S 





0 

B 


Here X and S denotes the actual values of the adjusted and unadjusted parameters, 
and 6X = X - X represents the difference between the actual and a priori values 

Qm Ct 

of the adjusted parameters. 


The normal matrix for all parameters is given by 


T 

] WG 


T 

A WA 

T 

A WB 

T 

B WA 

T 

B WB 


where the weight matrix W is the inverse of the observation noise covariance matrix 
and is assumed to be diagonal, i.e. 

= E(e£^0 = 

where E denotes expected values and is the ith observation noise sigma. 

The Bayesian least squares estimate of X is given by 

6X = (A'^WA+P“^)“^A\5y + (a\a+p“^)''^"^6X 
a ' a a a 

x^here 
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and 


= E[(X^-X)(X^-X)^1 

is the covariance matrix of the a priori estimate X^. . Substitutions yield, the 
error 

AX = 6X - SX = (A^WA+p“^)"^(a’^WB6S4-a\e-P“^6X ) 

3 - 3 3 

showing the three distinct components due to aliasing, measurement noise, and a 
priori uncertainties. 

Under the assumption that 6S, e, and 6X^ are uncorrelated errors, the covari~ 
ance matrix of the adjusted parameters is given by 

P = E[(X-2^X-X)'^] -= (A^WA+P”^)”^ 

a 

+ (a'^WA+p"^) (A^I'JB) (b'^WA) (a\a+p"^) 

3 ' o 3 

t 

where Vg is the diagonal covariance matrix of the unadjusted parameters, 

Pg = E(6S5S^) . 

The sensitivity matrix of adjusted to madjusted parameters is 

S ~ vrrjr = (A WA+P ) A WB 

dob a 

so that the aliasing error to the ith adjusted parameter due to the jth unadjusted 
parameter is 


The ith diagonal element of P, 
then given by 


a . . = 5 . . o„ . 
« « S. 


or the variance of the ith adjusted parameter is 


15 



a ^ = [(a\a+P ^)T^1 + y) 

XL a XX J ^ xj S, 

j ^ 

so that the standard deviation of the ith adjusted parameter is the root sum square 
(RSS) of the standard deviations due to data noise and a priori uncertainties 
and the standard deviations due to each unadjusted parameter. 
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5.0 MISSION STUDY CONFIGUMTION 


The ability of the "low- low" configuration to accurately recover geopotential 
fine structure of a given resolution (i.e. surface block size) depends on several 
factors : 

a. satellite altitudes 

b. satellite separation distances 

c. SST data types 

d. SST data rates 

e. SST data noise 

f. satellite state errors 

g. a priori values of the surface densities 

h. aliasing or orthogonality of the representation to the data type. 

Other studies [ (4), (9).] have investigated. the GRAVSAT SST missions for factors 
(a-e) with somewhat positive results.’ In an attempt to establish the feasibility 
of local fine structure recovery using "low-low" SST data when factors (f~h) are 
considered we have restricted the surface density blocks to be 5° x 5° and the 
polar satellite orbits to be circular at 250 km altitude. Moreover, our investi- 
gation and that of Schwarz and Hajela have shown that SST range data is much less 
sensitive to fine structure detail than SST range-rate data. All studies pre- 
sented' in Section 6.0 utilize SST range-rate data. 

The variation of the range-rate signal with satellite separation as the 
satellite pair passes over the centev Z8 density blocks of Figure 3 (the remaining 
blocks being neglected) is Illustrated in Figure 4-6 j where each block has been 
approximated by 16 mass points. The values used for the surface block densities 
in Figure 3 are those obtained by Schwarz after removing the full 12th order 
gravity field. Figures 4-6 clearly show the trade-off between increasing signal 
strength and decreasing resolution of detail as the satellite separation is 
increased. The aliasing of the signal signature over these same blocks by adja- 
cent blocks is illustrated in Figure 7, where now all 72 density blocks of Figure 
3 have been included in the 'force model. Comparison with Figure 5 shows the 
signal signature over the central region to have only slight alteration, while 
the signal toward the edges of the 36 central blocks is considerably altered. 

This indicates that on a single pass, the effect of unaccounted-for density 
blocks severely aliases nearby blocks, but that the effect rapidly falls off with 
increasing distance. The short arc error analysis studies of Section 6.0 clearly 
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demonstrate this phenomena. The situation fox long arcs is considerably different 
The presence of unaccounted-for blocks perturbs the satellite orbit and the effect 
builds up through the dynamics over long arcs so that distant blocks can cause 
serious aliasing. This effect can be seen in comparing Figures 7 and 5 by noting 
that the two curves have different values and slopes at the end of the pass. 

The long term perturbation of the orbit is clearly illustrated in Figure 8. 
Here the "low-low" configuration is in a circular equatorial two-body orbit at 
250 km altitude, and perturbed through four orbital revolutions by a single 1 mgal 
5° X 5° surface density block located on the equator. This shows that any 
unaccounted for block which perturbs the satellite orbit will alias the signal. 

The extent of the aliasing will be demonstrated in the long arc error analysis 
studies of Section 6.0. 

Our computer analysis using ERODYN and GEOMAP in the estimation mode has 
shown that the optimum satellite separation is approximately 6° , or 700 km for 
the 5° X 5° surface blocks at 250 km altitude. Consequently, all error analysis 
results presented in this report will assume this separation. Moreover, our 
ntimerical studies have shown that for error analysis purposes the numerical 
approximation for the surface density force model evaluation offered by GEOMAP 
Option (2) '(Section 3.0) with 9 mass points per block inside of a geocentric 
angle of 10° is of sufficient accuracy for the satellite configuration used and 
represented a substantial savings in computer time over Option (1) . 

The data rate assumed in this study is one observation every five seconds 
along the satellite path and four satellite passes per surface density block. 

This presents no problem for independent short arc solutions (less than one 
satellite period) since each arc normal matrix is generated separately and then 
concatenated, but for long arcs would require an excessive amount of computer run 
time to integrate the satellite orbits long enough to obtain the data coverage 
needed over a geographically localized set of density blocks. This problem was 
circumvented for the long arc- simulations by slowing the rotational rate of the 
earth by a factor of 18 so that the satellite ground track precessed approximately 
1.2° per revolution, producing the desired data coverage. 
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All data arcs consist of ascending satellite passes. This introduces some 
north-south asymmetry into the results, particularly for the short arcs, in that 
the. variational equations for the block parameters are all propagated from zero 
initial conditions at epoch at the south edge of the set of blocks under consid- 
eration. northward along the data arc. 

For final checkout of the GEO>IAP/ERODYN system, the GEOMAP program was 
operated in the estimation mode for long arc Tecovevies of the center 36 blocks 
of Figure 3. Data was generated using only a two-body force model and the 36 
perturbing density blocks. (This is, in fact, a global recovery in that no other 
blocks existed in the force model to alias the solution.) Recovery using perfect 
range— rate data (no noise) , zero a priori density values , 10 mgal a priori surface 
density sigmas and perfect knowledge of the satellite states yielded three digits 
of accuracy (approx. .01 mgal) in the estimated surface densities, giving confi- 
dence as to the correctness of the GEOMAP software. Adding random noise with 
sigmas of .005 cm/sec and .05 cm/sec to the data resulted in estimated surface 
densities with 2 digits (approx. 0.1 mgal) and one digit (approx. 1.0 mgal) of 
accuracy, respectively. Error analysis results using ERODYN established consis- 
tency with these errors of the GEOMAP least squares parameter estimation. 
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6.0 ERROR ANALYSIS STUDY 


From the example of aliasing in Section 5.0 it is clear that if a local set 
of density blocks is to be estimated the blocks near the edges will be badly 
corrupted by nearby unadjusted blocks, whereas blocks near the center of the set 
may be sufficiently far removed from the unadjusted blocks that their estimates 
may be acceptable. Thus, to estimate a given set of blocks it will be necessary 
to simultaneously estimate the set of interest plus a surrounding set of "buffer" 
blocks, whose estimates xjill be discarded due to aliasing. Moreover, the accuracy 
of the estimation of the center set of blocks of interest will depend on the 
region covered .by observations in the estimation procedure. To determine the 
relationships between estimation accuracy and the estimation region, data region, 
data noise, a priori uncertainties in adjusted and unadjusted density blocks and 
a priori imcertainties in the satellite states we utilize the method of covariance 
error analysis outlined in Section 4.0. It must be stressed that in a covariance 
analysis the least squares estimation is postulated and not actually performed; 
only the covariance matrix for the estimator is calculated. Parameters in the 
adjust or "solve for" mode are assumed to be estimated in the postulated least 
squares adjustment, while parameters in the unadjust or "consider" mode are 
assumed to affect the functionality between the adjusted parameters and the ob- 
servations but are not estimated in the postulated least squares adjustment. In 
several instances to be presented later the error analysis results somewhat strain 
intuition. However, it must be borne in mind that the error analysis problem 
treated here involves a large set of parameters, with many in the "consider" mode, 
and intuition can be misleading. 

Two basic sets of 5 ° x 5 ° equal angular density blocks are iised in the study: 

1. 225 blocks arranged in a 15 block by 15 block square centered on the 
equator with data coverage over the center 25 blocks. 

2. 289 blocks arranged in a 17 block by 17 block square centered on the 
equator with data coverage over the center 49 blocks. 

The arrays of blocks were set to these large dimensions so that the results ob- 
tained for subsets of blocks in the interior would not be compromised by possible 
non-physical "edge" effects. 
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The placement of the block arrays on the equator minimizes any "area effect" 
due to decreasing block areas with increasing latitude. The data coverage over 
the twenty-five blocks is accomplished with 21 ground track passes, while the 
49 block region required 31 passes. The data passes are symmetrically defined 
by the ground track of the midpoint between the two satellites, figures 9 and 
10 display the ground tracks over which data is taken for the two sets of blocks. 

The total parameter set considered for the error analysis consists of six 
state parameters for each satellite for each arc plus the entire set of surface 
density parameters- The sensitivity matrix C, consisting of the partial deriva- 
tives of the observations with respect to the total parameter set, is evaluated 
along a nominal two body trajectory hy setting the a priori values of all surface 
densities to zero. The resulting normal matrix is adequate for linear error 
analysis purposes and provides a substantial saving of computer time. 

The error analysis results were generated under the following assumptions: 

1. The two "low- low" satellites were configured in circular polar orbits 
at 250 km altitude with a 6° geocentric angle separation. 

2. Satellite-to-satellite range-rate data was generated at 5 second 
intervals (15 data points per 5° block) with a noise of .05 cm/sec. 

3. Data passes x^ere at approximately 1.2° spacing, giving approximately 
four data passes per block. 

4. All data was generated on ascending satellite passes. 

5- A priori satellite epoch state errors reflect the orbital knowledge . 
from other tracking means. No ground tracking data was included in 
the study. 

The error analysis results are organized into four sets: 

SET I . A single long data arc over the 225 block set. The twelve 

satellite state parameters are adjusted together with various 
subsets of blocks (Tables I.l through 1.6). 
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FIGURE 9 21 ascending short data arcs over 25 5“ x 5° density 

blocks. Ground track spacing is approximately 1.2°. 
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FIGURE 10 31 ascending short data arcs over 49 5° x 5*^ 

density blocks . Ground track spacing is 
approximately 1.2°. 
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SET II ♦ Twenty-one short data arcs over the 225 block set. All 252 

satellite state parameters are adjusted together with variotis 
subsets of blocks (Tables II. 1 through 11.20). 


SET III . Thirty-one short data arcs over the 289 block set. All 372 

satellite state parameters are adjusted together with various 
subsets of blocks (Tables III.l through III. 7). 

SET IV . Both a single long data arc and twenty-one short arcs over the 

225 block set. Aliasing is ignored and the satellite states are 
assumed perfectly known (Tables IV. 1 through IV. 4) . 


In the tables of solution sets I-IV displaying error sigmas, the dashed 
perimeter encloses the data region, while the solid perimeter encloses the set of 
adjusted density blocks. The two numbers within each adjusted block -represent 
the total RSS error sigma for that block over the error sigma for that block due 
to aliasing by unadjusted blocks in milligals (mgal) , i.e. 


'^TOTAL ^ ‘^ALIAS 


The error sigma due to" data noise and a priori parameter uncertainties is then 


'^NOISE 



2 

^TOTAL 


2 

'^ALIAS 


The numbers in unadjusted blocks represent the aliasing contribution of that block 

to the center adjusted block. The signs on these aliasing contributions reflect 

SAX 

the influence of the sensitivity matrix S = In the tables of solution sets 

I-IV displaying correlations, the number within each adjusted block is the correla- 
tion of that adjusted block with the center adjusted block. 


SET I; 

The purpose of this set of studies is to analyze the effect of aliasing on 
long arc solutions. Long data arcs for this problem have not been treated by 
other Investigators. As indicated in Section 5.0, to obtain reasonable computer 
run times, the rotational rate of the earth was slowed so that consecutive orbital 
ground tracks gave the proper spacing. This will certainly yield opti-mistic 
results in the error analysis since the dynamics will be correlated with a smaller 
number of unadjusted blocks than would be the case if the ground track passed 
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over larger areas of the globe. 


The normal matrix obtained for this study, as for all of Sets I-IV, dis- 
played poor conditioning for the recovery of the satellite epoch state parameters. 
Other studies have avoided this problem by not estimating them or assuming them 
to be perfectly knoxm (or, equivalently to have very small a priori sigmas) . 
However, we have found that when the satellite epoch state parameters are placed 
in the unadjust mode with sigmas as small as .001 meters/. 001 meters/sec for the 
positions and velocities, respectively (these numbers were chosen for convenience, 
the requirement being only that they be small) , their aliasing error completely 
swamps the estimation of the density blocks. The satellite epoch states, unless 
perfectly known, must therefore be estimated along with the density blocks. 

As could be guessed from the relative range-rate data type, the normal matrix 
is poorly conditioned for recovery of radial and cross track epoch state components, 
but is well conditioned for the recovery of along track epoch velocity components. 

In solution sets II and III this will be pursued by properly transforming sigma 
inputs in an HCL (radial, cross track, along track) coordinate system into a non- 
diagonal submatrix of the a priori covariance matrix. Otherwise all sigmas will 
be in cartesian coordinates in the diagonal a priori covariance matrix. In this 
solution set we set the epoch cartesian component sigmas to .OOJL_meters and .001 
raeters/sec for position and velocity to give strong diagonality to the matrix 

(A^WA -b P“^) 

to avoid the numerical difficulties in inversion. 

An estimation region which is exactly coincident with the data region (no 
buffer layer) is displayed in Table I.l. The a priori sigmas of the density 
blocks are 1 mgal. The recovery is seen to be very poor, with the a posteriori 
sigmas along the boimdary larger than the a priori value of 1 mgal. The errors 
are dominated by the aliasing contribution with the along track (N-S) being more 
severe than cross track (E-W) since the satellite orbits do not pass over the 
blocks outside of the longitude interval defined by the data region. An important 
point is the slowness with which the aliasing contribution decreases in the N-S 
direction. 
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The effect of adding a single "buffer layer" is sho^m in Table 1.2. The 
total error in the center of the estimated region decreased, but the aliasing 
error is still the dominant effect, showing a slow decrease with distance in the 
(N-S) direction. Table 1.3 displays the effect of adding two buffer layers. 

Only a slight improvement is obtained, with the (N-S) aliasing s till dominating 
the estimation accuracy. The correlations for the adjusted sets in Table I.l 
and 1.2 are displayed in Tables 1.4 and T. 5. The high correlations along track 
(N“S) reflect the fact that the relative range-rate data type is oriented along 
the satellite orbit. 

The implication of these results is that if there were blocks distributed 
over the globe, any block which the satellite pair passed over could alias the 
adjusted blocks. Table 1.6 clearly illustrates this fact. Here a band of 5°x 
5“ blocks completely encircles the globe passing through the poles. The data 
arc here consisted of a single long arc of 4 orbital revolutions with the earth 
slowed to give the proper data spacing. Again we see that the aliasing dominates 
.the error of adjusted blocks. As before we see the aliasing contribution of the 
unadjusted blocks falling off very slowly, then increasing to a maximum on the 
opposite side of the earth. 

This shows the low-low SST data type is unacceptable, due to severe aliasing, 
for recovery of localized density blocks for long arcs which pass“over unadjusted 
blocks. 


33 




LO 

4 ^' 


TABLE I.l 


■ 

lii 

r— — 





“.10 

B 

■ 

B 

■ 

■ 




1 




■ 


-.20 

1 

I 


SiS 




■ 

■ 






-.22 

B 


B 

B: 


• 


■ 

■ 






-.25 

B 

B 

B 

B 




■ 

■ 

■ 



H 

■ 

-.30 

B 

B 

B 

B 





■ 


■ 


I 

J 1 . 2 C 

.lit?. 

.97 


.92 

— 1 
JL^I 
1.24 1 


E 




In 

■ 

■ 

• 

■ 

i .99 

.02 

— 

olSL 

.58 

■' - • 

.60 

o 2 il | 
.95 1 

HI 

nii 




.02 

.02 

.02 

.001 

-.1 

m 

J12. 

.81 

.liss. 

.01 


uni 

1.03 I 

-.1 

.001 

.02 

,02 

.02 

■ 

■ 




* a 2 A 
j .94 

jjfii 

.GO 

, 5 G 

,62 

.97 j 

B 


• 








lizsa 

,. 9 S 

.93_ 

jjn 

*99 

Bl 

E 











■ 

-.30 


■MM 






■ 

• 





■ 

m 


B 

B 







■ 

■ 

■ 

■ 

■ 

B 


B 

B 

B 





■ 

■ 


■ 

■ 

■ 

B 


B 

B 

B 






li 


■ 

■ 

■ 

B 


B 

B 

B 


V* 



A Posteriori Density Sigmas, One Long Data Arc (Data Noise = .05 cm/eec) 
Adjusted parameters and a priori sigmas: 25 density blocks, a = 1 mgal 

Satellite States, *^pos~ *001 

°^VEL~ 

Unadjusted parameters and uncertainties; 200 density blocks* , o = 1 mgal 
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TABLE 1.2 A Posteriori Density Sigmas. One Long Data Arc (Data Noise = .05 cm/sec) 

Adjusted parameters and a priori sigmas; 49 density blocks, 0=1 tngal 

Satellite States, OpQg= .001 

Unadjusted parameters and uncertainties: 176 density blocks, 0=1 mgal 
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TABLE 1.3 A Posteriori Density Sigmas, One Long Data 
Adjusted parameters and a priori sigmas; 

.Unadinsted_jELa3"^nip1-r>rR and 11nf'pr^a•f nfi nq ; 


Arc (Data Noise = .05 cm/sec) 

81 density blocks,c 0 => 1 mgal 

Satellite .States, '^pos~ *001 

; a „ = .001 

VEL 
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TABLE 1.4 



Correlations of Adjusted Blocks with Center Block* One Long Data Arc 
(Data Noise - .05 cm/sec) 


Adjusted parameters and a priori sigmas; 25 density blocks, a 1 mgal 

Satellite States, o .001 

ir iJo 

°VEL" 

Unadjusted parameters and uncertainties; 200 density blocks, a = 1 mgal 
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■ TABLE 1,5 Correlations of Adjusted Blocks with 'Center Block,' One Long Data Arc 
(Data Noise - .05 cm/sec) 

■ Adjusted parameters and a priori sigmas: 49 density blocks, a = 1 mgal 

Satellite States, Op^g .001. 
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TABLE 1.6 


A Posteriori Density Sigmas One Long Data Arc (Polar Orbit over 4 
revolutions) 

Adjusted parameters and a priori sigmas: 6 density blocks, a = 1 

Satellite States c^Qg = «■< 

NOTE: Aliasing errors are with respect to adjusted block marked by pointer 


















































SET II: 


Solution Set II consists of 21 independent short data arcs over the 25 center 
blocks. All results display the' N-S asymmetry discussed in Section 5.0 due to 
the use of only ascending arcs with epochs along the south edge of the data region. 
This results in shifting the smallest error sigmas from the center of the adjusted 
set of blocks toward the south boundary in several cases. However, the center 
block remains the Ttiaiu focus of attention in these comparative studies for quanti- 
tative aliasing contributions. 

1 

Tables II.1-II.3 are to be compared directly with I. 1-1.3. In both sets, the 
a priori error sigmas on the density blocks and the satellite states were 1 mgal 
and .001 m/.OOl in/s (virtually no error, .001 chosen for convenience), respectively. 
The important point to note is that the aliasing error is no longer the dominat- 
ing error source, as the aliasing contribution of unadjusted blocks decreases 
rapidly with distance from the adjusted set. Moreover, the additive "buffer 
layers" of Tables II. 2 and II. 3 show dramatic decreases in the aliasing error 
contribution to the center block. The dominant error source is now the "noise 
plus a priori" contribution arising from 

•p 

(A WA + ) . 

The complicated dependencies of the components of total error on the a 'pvioT-i 
sigmas for the satellite epoch state and the density blocks are investigated in 
Tables II. 5 through 11.20. 

Tables II. 5, II. 6 and II. 7 where the' a priori errors on the satellite epoch 
states are systematically increased are designed to be compared with Table II. 1. 
Table II. 5 displays the surprising result that an increase in the satellite epoch 
state a priori to a = .002 (meters and meters/sec) resi^lts in a decrec^e in both 
the total error and the aliasing contribution to the center adjusted block. It 
must be noted, however, that there are increases in error to some of the adjusted 
blocks and that the total error to aVt adjusted parameters (satellite epoch states 
plus the 25 adjusted density blocks) increases. Table II. 6, with an increase in 
the a priori epoch sigmas to a = .01 shows an inovease in the total error sigroa 
of the center block (but a decrease in the alias error) while some of the blocks 
in the northern sector again show a decrease. As indicated previously, the 
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covariance matrix is poorly conditioned for the estimation of the radial and 
cross track components of the satellite states, as would be expected for the SST 
range-rate data. The inversion of the normal matrix plus the inverse a priori 
covariance matrix broke do^m due to loss of significant digits when the a priori 
velocity sigmas were increased beyond a = .01 m/s for the Cartesian components. 

In an attempt to utilize realistic sigmas, the option' to input satellite epoch 
sigmas in a radial, cross track, along track (denoted HCL) coordinate system was 
implemented. The maximum satellite state error sigmas permissible before inver- 
sion failure (with triple precision arithmetic for the inversion algorithm) was 
detetnoined to be 10 meters in the position components and .01 meters/sec in the 
radial and cross track velocity components. As would be expected for the data 
type, the along track velocity showed a high recoverability, with error sigmas 
up to 100 meters/sec showing no Inversion difficulties and the a posteriori estimated 
error sigma less than 1 meter/sec. Tables II. 7 and II. 8 show the effect on the 
adjusted blocks of increasing the along track a priori velocity sigma values to 
2 meters/sec and 5 meters/sec, respectively, with 10 meters in position and .01 
meters/sec in the radial and cross track velocities. 

The effect of the a priori sigmas of the density blocks is investigated in 
Tables II.9-II.il where the density block sigmas are set to 50 mgal and the esti- 
mation configurations of Table II. 1 through II. 3 are repeated. This set of solu- 
tions clearly demonstrates the importance of adding buffer layers. Wien the 
adjusted region covers the data region and one buffer layer aromd 
the data region the aliasing due to the large uncertainties in the unadjusted blocks 
is the dominating effect. The addition of two buffer layers in Table II. 11, 
however, drastically reduces both the total error of the interior blocks and the . 
fraction of the total error due to aliasing. It should be noted in comparing 
Tables II. 1 and II. 9 that although the estimation strategy of adjusting all blocks 
in the data region with no buffer layers is a poor one, it does show some improve- 
ment with an assumed a priori knowledge of all density blocks to 1 mgal. The 
assumption of a 50 mgal a priori uncertainty in the density blocks, however, re- 
sults in a worsened solution. 

The effect of systematically increasing the a priori surface density sigmas 
for a one layer buffer zone and satellite epoch state sigmas of 10 meters in 
position, .01 meters/sec in radial and cross track velocities, and 2 meters/sec 
in along track velocity is presented in Tables 11.12 through 11.17. The case of 
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100 mgal a priori uncertaiaty of Table 11.15 shows the aliasing error again 
becoming the dominant effect. Table 11.16 shows the effect of having good a 
priori knowledge of the’ adjusted blocks (a =1 mgal) and little a priori 
knowledge of the unadjusted blocks (a = 50 mgal). The recovery is destroyed by 
the dominating alias error. The opposite extreme of poor a priori knowledge of 
the adjusted set (50 mgal) and good a priori knowledge of the unadjusted set 
(a = 1 mgal) is shown in Table 11.17. As would be suspected, this recovery 
suffers little from aliasing errors. The case for two buffer layers and a priori 
surface density sigmas of 50 mgal is shown in Table 11.18. Comparing this result 
with the single buffer layer of Table 11-14 shows a much less dramatic improve- 
ment in estimated accuracy for the adjusted blocks than the similar cases with 
small a priori sigmas on the satellite states (Tables II. 10 and II. 11), indicating 
the importance of some amoimt of a priori knowledge. 

The effect of reduced data noise is displayed in Tables 11,19 and 11.20 
where the noise sigma is set to .01 cm/sec, while the satellite state sigmas are 
kept small (.001) and the surface density block sigmas are set at 1 mgal. These 
results are to be compared with Tables II. 2 and II. 3 where the data noise is .05 
cm/sec. The total error decreased in both cases 11.19 and 11.20, while the 
aliasing contribution increased slightly. The dominating error source in these 
cases is the a priori imcertainties in the adjusted parameters. 
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TABLE 11,2 A Posteriori Density Sigmas • ' ' . 

21 Short Data Arcs (Date Noise - ,05 cm/sec) 

Adjusted Parameter and a priori sigmas; 49 density blocks, a = 1 mgal 
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TABLE II. 7 A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise = .05 cm/sec) 

Adjusted parameters and a priori sigmas; 25 density blocks, a = 1 ragal 

Satellite States, Op^g = 10 
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TABLE II. 8 A Posteriori Density Sigmas ' • ’ ■ ■ 

- 21 Short Data Arcs (Data Noise = .05 cm/sec) . , 
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TABLE 11.12 ' A Posteriori' Density Sigmas • , ‘ 

21 Short Data Arcs (Data Noise = ,05 cm/sec) 

Adjusted' parameters and a priori sigmas; 49 density blocks, a = 1 mgal 

Satellite States, o =10 
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A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise = .05 cm/sec) 
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TABLE 11.14 A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise = .05 cm/sec) • 

Adiusted parameters and a priori sigmas; 49 density blocks, . = 50 meal 

Satellite States, cr^Qg = 10 

a^gj^EHCL] = (.01, . 

Unadjusted parameters and uncertainties; 176 density blocks, cr =• 50 mgal 
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TABLE 11.15 A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise = .05 cm/sec) 

Adjusted parameters and a priori sigmas: 49 density blocks, a - 100 mgal 

Satellite States ^fpQg “ 

o^^EHCL] = (.01, . 

Unadjusted parameters and uncertainties: 176 density' blocks , o = 100 mgal 
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TABLE 11.16 A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise = .05 cm/sec) 

Adjusted parameters and a priori sigmas: 49 density blocks, o == 1 mgal 

Satellite States, Op^g = 10 

Unadjusted parameters and uncertainties: 176 density blocks, o = 50 mgal 
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TABLE 11.17 A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise - .05 cm/sec) 

■ Adjusted parameters and a priori sigmas.: 49 density blocks, o = 50 mgal 

Satellite States, 10 

JrUo 

. 0^gj^[HCL] = (.01, .01, 2) 

Unadjusted parameters and uncertainties: 176 density blocks. O' == 1 mgal 
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TABLE 11.19 A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise = .01 cm/sec) 

Adjusted parameters and a priori sigmas: 49 density blocks, a = 1 mgal 

Satellite States, cfp^g = .001 

cf = .001 

VEL 

Unadjusted parameters and uncertainties: 176 density blocks, o = 1 mgal 
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TABLE 11.20 


A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise - .01 cm/sec) 

Adjusted parameters and a priori sigmas: 81 density blocks. 


Satellite States, 0^ 


Unadjusted parameters and uncertainties: 144 density blocks, a 





















































































SEX III. 


Solution set III utilizes a 17 x 17 rectangular array of 5 ° x 5 ° density 
blocks with 31 independent data arcs giving data coverage over the center 49 
blocks. The adjusted parameters in. this set consist of the 372 satellite epoch 
state parameters plus various configurations of adjusted blocks over the data 
region. The principal design of this solution set is to display the effect of 
increasing the data region relative to that of Solution set II for both small 
(1 mgal) and large (50 mgal) a priori density block sigmas. All solutions have 
small a priori Cartesian error components on the arc states (.001 meters and 
.001 meters/sec for position and velocity, respectively). 

Table TII.l displays the disastrous effect of estimating a subset of blocks 
within the data region, as most of the a posteriori sigmas of the recovered blocks 
are greater than their a priori values. It is clear that processing data which 
is over the unadjusted blocks and therefore strongly influenced by them should 
result in more significant aliasing errors on the adjusted set of blocks. 

Tables III. 2 through III. 4 display estimation regions with no buffer layer, 
one buffer layer and two buffer layers, respectively, for the case of one mgal a 
priori error sigma on the density blocks. The recovered sigmas of the central 
adjusted region compare very closely with the same configurations of Solution Set 

II (Tables II. 1 through II. 3) with a smaller region of data coverage. Tables Til. 5 
through III. 7 repeat the same estimation region configurations, but with a large 

a priori uncertainty of 50 mgal on the density blocks. -Comparing with the similar 
sets in Solution Set II (Tables II. 9 through II. 11) shows that increasing the 
data region with a poor estimation strategy (no buffer layer or only one buffer 
layer) and a large a priori uncertainty on the density blocks can result in con- 
siderably increased aliasing errors. The region covered by data in Solution Set 

III is nearly twice that of Set II, so that there are correspondingly more 
unadjusted blocks to contaminate the estimation- accuracy when a poor strategy is 
used. It should be noted that Table III. 7, with two buffer layers, again compares 
very closely with Table II. 11, which has the same configuration but a smaller data 
region. 
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TABLE m.l A Posteriori Density Sigmas 

31 Short Data Arcs (Data Noise = .05 cm/sec) 

Adjusted parameters and a priori sigmas; 25 density blocks, o = 1 mgal 

Satellite States, = .001 

'"VEI, - 

Unadjusted parameters and uncertainties: 264 density blocks, a = 1 mgal 































o> . 

Ln 



■ 

i 

i 

i 

■ 

■ 

i 

■ 

.0045 

i 


■ 

■ 


■ 

■ 

i 

■ 

■ 

■ 

■ 

■ 

i 

■ 

■ 

.0076 

■ 



■ 

■ 



■ 

■ 

■ 

■ 

■ 

■ 

■ 


■ 

.0131 



■ 

i 



i 

■ 

■ 

■ 




■ 

■ 

■ 


j/h 



■ 

1 


1 

■ 

■ 

n 

■ 

[I 

■ 

■ 

■ 

■ 

jjj^Q 

H 


■ 

■ 

■ 


■ 

III 

■ 

■ 

■ 

■ 

■ 

im 

.755 

.679 

■1.653. 

.574 

,666 

.722 

jugnii 

■ 

i 




m 

.889 

'.687 

.598 

.564 

.651 

.840 

■ 

■ 

■ 

■ 

■ 

.993 

,750 

.668. 

.559 

.650 

.657 

.733 

.94 5 


1 

■ 

1 

■ 

.933 

,658 

.536 

.546 

.640 

.880 






1-106 

,712 

,618 

BH 

■ 614. 

,688 

,943 

■ 


■ 

■ 

■ 


.573 . 

.444 


.439 

.543 

.852 


Bl 


iQm 

B 

I . 1 ?4 
1.053 

.526 

■t.566. 

,351 

1^9 

,553 

-■-6.18. 

.445 

,804 

s 


.0007 



.345 

.340 

.783 




1 

■ 


,661 

.530 

,504 

,506 

.402 

.784 

■ 


■ 


■ 

t-ul 

.530 

,342 

.300 

.305 

.686 





■ 

..afti’K 

,664 


.472 

HH 

.561 

..703 

■ 


■ 

■ 

■ 

.790 

.583 ■ 

.336 

.361 

.460 

.632 





! 

■ ■?94 

WSSk 


,472 

9SS 

.548 

.,708 




1 

■ 


.552 

• 37^ 

r3£ij 

.471 

•6S1 






■ 

■ 

11 


■ 

■ 

m 

■ 


■ 

1 

■ 






■ 

■ 

■ 


■ 

i 

m 




i 

■ 


■ 




i 


■ 


i 

H 

■ 

■ 



i 


■ 

1 




■ 


■ 


i 

■ 

■ 

■ 

m 

■ 


■ 

i 

1 





1 

1 


■ 

■ 

■ 

1 

■ 

1 

■ 

H 


TABLE III. 2 A Posteriori Density. Sigmas 

31 Short Data Arcs (Data Noise “ .05 cm/sec) 

Adjusted parameters and a priori sigmas: 49 density blocks, 0=1 mgal 

Satellite States, cTp^g = .001 


Unadjusted parameters and uncertainties: 240 density blocks, o "-1 mgal 
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TABLE III. 3 


A Posteriori Density Sigmas 

31 Short Data Arcs (Data Noise '= .05 cm/sec) 

Adjusted parameters 'and a priori sigmas: 81 density blocks, a “ 1 mgal 
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A Posteriori Density Sigmas 

31 Short Data Arcs (Data Noise = .05cm/sec) 
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A Posteriori, Density Sigmas 

31 Short Data Arcs (Data Noise - .05 cm/sec) 

Adjusted paraijieters and a priori sigmas: 121 density blocks, a = 50 mgal 

Satellite States ~ -001 


‘^VEL 


.001 


Unad j us ted parameters and uncertainties 
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SET IV. 


Solution Set IV displays the error analysis for an estimator which ignores 
aliasing effects and assumes the satellite states to be perfectly known. This in 
some measure serves as a bridge between the results of Schwarz, Hajela and 
Solution Sets I and II, showing clearly how misleadingly optimistic the results 
of error analysis studies can be when important aliasing errors are neglected. 

Table IV. 1 examines the long data arc of Solution Set I with data over the ■ 
center 25 block region. The twelve satellite state parameters are put into the 
unadjust category along with all density blocks outside of the data region .with 
a zero a priori sigma. The total error is then the "noise only" contribution; 
the aliasing error is zero. The two numbers inside the twenty-five adjusted 
density blocks represent the recovered sigmas for the two cases of 1 mgal and 
50 mgal a priori sigma for the adjusted blocks, respectively, i.e. 

c 

*^1 mgal' ^ ^50 mgal 

Unlike the cases of Sets I, II and III, the lack of a priori knowledge of the 
adjusted blocks has little impact on the estimation accuracy. Table IV. 2 dis- 
plays the correlation coefficients of the adjusted blocks with the center block' 
for the case depicted in Table IV. 1. The N-S and E-W patterns closely resemble 
the results of Schwarz. 

Tables IV. 3 and IV. 4 correspond to the error analysis of Tables IV. 1 and IV. 2 
with the exception that 21 independent data arcs are used in place of one long 
arc. Generally the same patterns result with the exception of the asymmetry 
introduced in the southern-most row -due, as previously discusssed, to the varia- 
tional equations for the density block parameters all starting from zero initial 
conditions at epoch. It should be noted that the long arc solution of Table IV. 1 
displays a superior estimation accuracy to the short arc solution of Table IV. 2, 
even though the data coverage is the same. The reason for this, of course, is 
that -the short arcs are assumed independent, while the consecutive passes of the 
long arc are correlated through the orbital dynamics. 

The solution displayed in Table IV. 3 corresponds quite closely with the 
solution 3.1 of Schx<rarz for 5“ x 5° blocks with the IoX’T-Iow configuration at 300 km 
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altitude when proper corrections are made for the different values used for the 
data noise and for the differences in data rates and data spacing. Noting that 
Schwarz used approximately one data arc per 5® x 5 ° block, an observation approx- 
imately every 30 seconds and a data noise of .005 cm/sec, a factor to relate his 
estimated density parameters to those of Table IV. 3 would be 



Thus his estimated uncertainty of approximately .2 mgal should correspond to 
approximately .4 mgal for our results of IV. 3. As the satellite altitude used 
in Table IV. 3 was 250 km compared to 300 km used to Schwarz in solution 3.1, we 
would expect his values to be slightly larger than ours. 
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TABLE IV. 1 A Posteriori Density Sigmas 

One Long Data Arc (Data Noise = .05 cm/ sec) 

Adjusted parameters and a priori sigmas: 25 density blocks, 0 = 1 mgal 

50 mgal 

Unadjusted parameters: None 



























■ 

■ 
















■ 













■ 


■ 



■ 

■ 

IHi 








■ 



• 

■ 

1' 

■ 

■ 








■ 


■ 

■ 

lilH 

■ 

■ 

■ 










,-.0l 

.03 

-.04 

.09 

1 

-.07 j 

III' 








■ 


l-.os 

.13 

-.31 

. 13 
.14 

-.05 , 

... 1 

■ 








■ 


1 uM 

1 *27 

-.61 

1 

-.60 

-.61 

^ 1 
.27 1 

1 








■ 


I-. 05 
|-,05 

■H 


oil 

.13 

-.05 1 
-.05 1 





i 






pi.Oi 

j-. 0''* 

.08 

-.03 

.02 

1 

^.,^1 

1 




1 




■ 


H 

in 




■ 








■ 


■ 

■ 




■ 








11 

. 


■ 




■ 








■ 



■ 




■ 








■ 


H 





■ 




. • 


TABLE IV. 2 ■ Correlations, of Adjusted Blocks with. Center Block ■ . ' ' 

One Long Data Arc (Data Noise = '.05 cm/sec) 

Adjusted parameters and a priori sigmas: 25 density blocks., o = 1 mgal 

50 mgal 

Unadjusted parameters; None 
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A Posteriori Density Sigmas 

21 Short Data Arcs (Data Noise = .05 cm/sec) 

Adjusted parameters and a priori sigmas: 25 density blocks, a = 1 mgal 

50 mgal 


Unadjusted parameters: None 
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TABLE IV. 4 Correlations of Adjusted blocks with Center Block 
21 Short Data Arcs (Data Noise = .05 cm/ sec) 

Adjusted parameters and a priori sigmas: 25 density blocks, a = 1 mgal 

50 mgal 

Unadjusted parameters: None 



















7.0 SUJmRY 


This study applied covariance error analysis techniques to investigate 
estimation strategies for the low-low-SST mission for accurate local recovery of 
gravitational fine structure, considering the effects of aliasing. A 5° x 5“- 
surface density block representation of the high order geopotential was utilized 
with the drag-free "low-low" GRAVSAT configuration in a circular polar orbit at 
250 km altitude. 

Satellite-to-satellite relative range-rate data was used in the error analysis 
investigation, as the range data was found to be much less sensitive to fine 
structure recovery. A data noise of .05 cm/sec was assumed for most of the error 
analysis studies. No ground tracking data was included and satellite state 
accuracy requirements for the estimation of the density blocks were investigated 
in the error analysis by a priori error sigmas for the epoch states, both in 
Cartesian and HCL coordinates. The sensitivity of the data type to even very 
small a priori errors in the satellite epoch state velocities (.001 meters/sec) 
led to severe aliasing errors in the estimation of local sets of density blocks, 
forcing the satellite epoch state parameters to be included with the surface 
density parameters in the estimation process. 

The SST relative range-rate data type exhibits strong recovery ability only 
for' the along track components of velocity the normal matrix being poorly -con- 
ditioned for the recovery of radial and cross track components. This requires 
satellite orbital accuracies from other tracking means to a minimum of 10 meters 
in. position and .01 meters/sec in radial and cross track velocity components. 

The recovery of local sets of density blocks from long data arcs proved not 
to be feasible due to strong aliasing effects from non-estimated density blocks 
all along the ground track - even those very far removed from the set of estimated 
blocks. The error analysis for the recovery of local sets of density blocks using 
independent short data arcs demonstrated that the estimation strategy of simultane- 
ously estimating a local set of density blocks covered by data and two "buffer 
layers" of density blocks not covered by data resulted in almost negligible alias- 
ing errors on the estimation accuracy for blocks near the center of the recovered 
set due to unadjusted surface density blocks. 
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With an orbital accuracy of 10 meters in position, .01 meters/sec in radial 
and cross track velocity and 2 meters/sec in along track velocity obtained by 
other means of tracking, the low-low SST range- rate data, with a noise of -05 
cm/sec, will recover 5° x 5° surface density blocks to approximately 5 mgal 
and 17 mgal, respectively, when an a priori knowledge of the blocks to 15 mgal 
and 50 mgal is assumed. [Note: These values wheii multiplied by 2ir give the 

approximation to the corresponding values of gravity anomaly blocks.] This is 
to be contrasted to results (Solution Set IV) predicting recovery to approximately 
.2 mgal (40 cm geoid) when the error analysis does not consider aliasing errors 
from unadjusted density blocks and assumes the satellite states to be perfectly 
knoTm. 

We would conclude that the estimation of gravitational fine structure by 
recovery of local sets of density blocks utilizing SST relative range rate data 
x-;ith the "low-loxy drag-free satellite configuration is not a feasible approach, 
principally due to the very stringent orbital accuracies required x^ith many short 
data arcs. The low altitude satellite orbits would make the attainment of such 
accuracy from ground tracking extremely difficult. 

The possibility of supplying such tracking from a third high altitude satel- 
lite is a concept (the ’’high-low-lox/' configuration) which des erve s investigation. 


An investigation of the feasibility of a global recovery for the gravitational 
fine structure parameters, a concept summarily ruled out of hand by many authors 
due to the magnitude of the problem, also deserves attention. Such a procedure 
would eliminate the aliasing problem inherent with the estimation of local sets 
of parameters and relax the difficulty of obtaining dense coverage of very accurate 
ground tracking. 



78 



8 . 0 ACKNOX^LEDGEMENT 


It is with pleasure that the authors acknowledge the support of Dr. P. 
Argentiero of GSFC, whose error analysis techniques developed in his study of 
the recovery of local sets of gravity anomalies were used in this investigation. 
The assistance provided by Mr, Daniel Chin of BTS in modifying the ERODYN and 
GEODYH programs and by Mr. Tom Englar and Mr. Andy Jazwinski of BTS for interest- 
ing discussions is also acknowledged. 


'ORIGINM] PAGS IS 
GS POOR 'gUAXOT 


79 



9.0 IffiW TECHNOLOGY 


The effort under this contract consisted of the application of covariance 
error analysis techniques for investigation of estimation strategies for the 
"low-low" SST mission for accurate local recovery of gravitational fine structure, 
considering the effects of aliasing errors. Frequent reviews and a final survey 
for new technology were performed. It is believed that the mathematical and 
programming techniques and algorithms developed do not represent "reportable 
items", or patentable items, within the meaning of the New Technology Clause. 

Our reviews and final survey found no other items which could be considered 
reportable items under the New Technology Clause. 
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